Jet Quenching and Azimuthal Anisotropy of Large px Spectra in Non-central 
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Parton energy loss inside a dense medium leads to the suppression of large pr hadrons and can 
also cause azimuthal anisotropy of hadron spectra at large transverse momentum in non-central 
high-energy heavy-ion collisions. Such azimuthal anisotropy is studied qualitatively in a parton 
model for heavy-ion collisions at the RHIC energies. The coefficient V2{pt) of the elliptic anisotropy 
at large pr is found to be very sensitive to parton energy loss, ft decreases slowly with pr contrary to 
its low pt behavior where V2 increases very rapidly with pr- The turning point signals the onset of 
contributions of hard processes and the magnitude of parton energy loss. The centrality dependence 
of V2{pt) is shown to be sensitive to both size and density dependence of the parton energy loss 
and the later can also be studied via variation of the colliding energy. The anisotropy coefficient 
V2/e normalized by the spatial ellipticity e is found to decrease significantly toward semi-peripheral 
collisions, differing from the hydrodynamic results for low pr hadrons. Constrained by the existing 
WA98 experimental data at the SPS energy on parton energy loss, both hadron spectra suppression 
and azimuthal anisotropy at high pr are predicted to vanish for b >7-8 fm in Au + Au collisions at 
^/s =130-200 GeV when the hadron rapidity density per unit area of the initial overlapped region 
is less than what is achieved in the central Pb + Pb collisions at the SPS energy. 
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I. INTRODUCTION 

In the study of the dense matter that is created in 
high-energy heavy-ion coUisions, one crucial issue is the 
degree of thermaUzation through secondary scatterings. 
If thermaUzation has been, at least partially, achieved, 
there should be collective effects developed through the 
evolution of the system. Azimuthal anisotropy in hadron 
spectra or elliptic flow has been proposed as a good signa- 
ture of collective transverse flow in relativistic heavy-ion 
collisions 0|. Such an elliptic flow has been observed in 
experiments at the BNL/AGS |], CERN/SPS [§ and 
most recently BNL/RHIC by the STAR experiment [|) . 
The strength of the elliptic flow, V2, experimentally de- 
fined as the second coefficient in the Fourier decomposi- 
tion of the particle azimuthal distribution p] with respect 
to the reaction plane, was found jq] to be between the 
limits of low-density rescattering M and hydrodynamic 
expansion |l|,p|-pi)[, indicating approach to a higher de- 
gree of thermalization with increasing colliding energies. 
At RHIC energies, one expects to see v^ becoming closer 
to the hydrodynamic limit as the initial energy density in- 
creases and the lifetime of the initial dense matter is get- 
ting longer. Because system evolution will eliminate the 
geometrical anisotropy which generates the anisotropy in 
momentum space, the elliptic flow was argued |1 l|,p| to 
be sensitive to the early dynamics of the system. 

In hydrodynamic models, the strength of the differen- 
tial elliptic flow, V2{pt), increases almost linearly with px 
|0,|| (at Pt ~ 0, however, W2 increases quadratically with 
Pt) because of collective expansion. At large transverse 
momentum, the hydrodynamic model will likely cease to 



be a valid mechanism for particle production in high en- 
ergy nuclear collisions. Instead, particle production at 
Pt > 2 GeV/c will be dominated by hard or semihard 
processes. The hadron spectra in this region of phase 
space typically exhibit a power-law behavior and depend 
on how an energetic parton propagates through the dense 
medium created in the heavy-ion collisions. The corre- 
sponding azimuthal anisotropy of hadron spectra at high 
Pt would also have completely different behavior from 
hydrodynamic model results. 

Recent theoretical studies of fast partons propagating 
inside a dense medium all suggest a large energy loss 
caused by multiple scattering and induced gluon radia- 
tion |l4-[l4| . Of particular interest is the quadratic de- 
pendence of the total energy loss on the distance of prop- 
agation [|l3[ due to the non-Abelian nature of gluon ra- 
diation in QCD. Such a parton energy loss will cause the 
suppression of large pt hadron spectra |la,n6[ in heavy- 
ion collisions if the lifetime of the dense medium is long 
enough to influence the propagation of fast partons. The 
suppression will depend on the average distance that par- 
tons propagate inside the medium during the lifetime of 
the dense matter. Since the transverse distance depends 
on the azimuthal direction of the parton propagation in 
non-central heavy-ion collisions, one should expect the 
hadron suppression caused by parton energy loss to de- 
pend on the azimuthal angle with respect to the reaction 
plane, thus leading to azimuthal anisotropy in high pT 
hadron spectra. 

When the momentum transfer ^ for each of the few 
scatterings suffered by an energetic parton in a finite sys- 
tem is small so that h/pt <C 1, the effect of the elastic 



scattering on azimuthal anisotropy in large pr hadron 
spectra is negligible. Only parton energy loss can cause 
sizable azimuthal anisotropy for large pr hadrons. There- 
fore, azimuthal anisotropy and hadron spectra suppres- 
sion at large pt should accompany each other in high 
energy heavy-ion collisions. For a parton energy loss 
dE/dx that has a weak dependence on the parton en- 
ergy, the hadron spectra suppression and the accompa- 
nying azimuthal anisotropy should all decrease with px 
as we shall show in this paper. This is in sharp con- 
trast with the effect of hadronic interactions whose cross 
sections, especially for inelastic processes, increase with 
energy and thus lead to increased hadron suppression and 
azimuthal anisotropy at large px- Therefore, the pT de- 
pendence of hadron spectra suppression and azimuthal 
anisotropy is a unique signal of parton scattering and 
energy loss in the early stage of heavy-ion collisions. 

Azimuthal anisotropy due to jet quenching has been 
studied J|l|,|l3 before with the HIJING Monte Carlo 
model |19f| T^ In this paper, we will study the azimuthal 
anisotropy in hadron spectra at large px using a parton 
model in which one incorporates the parton energy loss 
via modified parton fragmentation functions [16 20 1 . We 
will study the sensitivity of the azimuthal anisotropy, its 
Pt and centrality dependence to the parton energy loss 
dE/dx and in particular the distance and density depen- 
dence of the parton energy loss. 

We should emphasize here that our study in this paper 
is only qualitative. First of all, the hard sphere geome- 
try we use does not give an accurate description of the 
spatial anisotropy e (or ellipticity). However, the ratio 
W2/e will reduce the sensitivity to models of geometry. 
Secondly, our theoretical understanding of parton energy 
loss is only qualitative so far. There are many uncer- 
tainties in the estimate of dE/dx in QCD. The purpose 
of our study is to demonstrate the effect of jet quench- 
ing in hadron spectra and the azimuthal anisotropy for 
a given averaged dE/dx. Furthermore, in our model we 
will neglect both longitudinal and transverse expansion. 
The longitudinal expansion will change the particle den- 
sity and thus the effective total parton energy loss ||2^] 
in the medium while the transverse expansion will influ- 
ence the spatial ellipticity of the medium as the system 
evolves with time. Since the parton energy loss in a QGP 
is estimated to be much larger than that in a hadronic 
gas llJ , we assume that the parton energy loss happens 
mostly in the partonic medium. Then an energetic par- 
ton is likely to travel to the outside of the partonic matter 
before the transverse expansion changes the spatial ellip- 
ticity of the partonic medium significantly. In this case, 
the expansion will just change the overall average parton 
energy loss. For a qualitative study in this paper, we are 
only concerned with the average parton energy loss and 
we will just consider the dependence of dE/dx on the 
average particle density in the initial overlapped region. 



II. HADRON SPECTRA AT LARGE Pt 

Hadron spectra in pp, pA and AA collisions have been 
systematically studied in a parton model m& . This model 
extends the collinear factorized parton model to include 
intrinsic transverse momentum and its broadening due to 
multiple scattering in nuclear matter. The value of the 
intrinsic transverse momentum and its nuclear broaden- 
ing are adjusted once and the model can reproduce most 
of the experimental data in pp and pA collisions |16| . In 
AA collisions, we model the effect of parton energy loss by 
the modification of parton fragmentation functions |2C| ] 
in which the distribution of leading hadrons from parton 
fragmentation is suppressed due to parton energy loss. In 
non-central collisions, the single inclusive hadron spectra 
in this parton model are given by. 



dypTdprdcj) ^ 

abed 



dHd'^rtA{r)tB{\h~r\ 



/ dxadxtd^kaxd^ hrgAikaT , Q'^,r)gB{kbT, Q^ 



Ib-r| 



fa/AiXa,Q^,r)f,,/BiXb, 
Dh/c{Zc,Q^,L{<j)))d<T 



dt 



|b-r|) 
{ab -^ cd). 



(1) 



The factor iir = 1.5 — 2.0 is used to account for higher 
order corrections. The nuclear thickness function tA{b) 
is normalized to J d'^btA{b) = A using the Woods-Saxon 
form of nuclear distributions. To take into account of the 
intrinsic transverse momentum and its nuclear broaden- 
ing, gA{kT,Q^,r) is assumed to have a Gaussian form 
and its width is parametrized to fit the pp and pA 
data. The impact-parameter dependence oi gA{kT, Q^,r) 
comes from the nuclear broadening of the intrinsic trans- 
verse momentum. The parton distributions in nuclei 
fa/Ai^a, Q^, t) are assumed to be factorized into the par- 
ton distributions in a nucleon and the impact-parameter 
dependent nuclear modification factor. 
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We will use the parameterization of the nuclear modifi- 
cation in HIJING [H. Constrained by the existing pA 
data, nuclear modification of parton distributions and 
Pt broadening produce about 10-30% increase in the pt 
spectra in central Au + Au collisions relative to pp [|6| . 
However, they will not give any azimuthal anisotropy in 
the hadron spectra. 

Eq. ffl) also applies to pp{p) collisions in which one sim- 
ply sets A = 1 without intrinsic pt broadening and nu- 
clear modification of the parton distributions. The par- 
ton fragmentation functions Dj^ii,{zc, Q^) are then given 
by the parameterization of e~^e~ data pl[ . As an ex- 
ample, we show in Fig. |] the comparison of the parton 



model calculation (solid line) with the experimental data 
H for pp collisions at ^i = 200 GeV. We used aK = 1.5 
factor in the calculation. The model describes the data 
very well down to pr ^ 1-5 GeV/c. Below this scale we 
expect such a parton model calculation to fail and non- 
perturbative physics to dominate. As an illustration, we 
fit the experimental data below pt=1 GeV/c with an 
exponential distribution (thin solid line), 
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with Tq = 0.2 GeV. The normalization of the fit is de- 
termined by drich/dr] = 2.3 in the central rapidity region 
and (7in = 42 mb at y/s = 200 GeV. One can notice 
that the data differ significantly from the exponential fit 
already at around pr — 1-5 GeV/c. Beyond this scale, 
the spectra have a power-law behavior characteristic of 
hard parton scattering. We can therefore use the parton 
model to describe the hadron spectra for pT > 2 GeV/c. 
In Fig. ^ we also plot the hadron spectra for pp collisions 
at -^/s = 130 GeV. The spectrum differs very little from 
that at v^ ~ 200 GeV at intermediate pt and is only 
about a factor 2 lower at pt = 7 GeV/c. 
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FIG. 1. Transverse momentum spectra of charged hadrons 
in pp collisions at y/s = 200 (solid) and 130 GeV (dot-dashed) 
as compared to UAl data |23] on pp collisions at ^/s = 200 
GeV. The thin solid line is a parameterization of the data 
below Pt < 1 GeV/c. 

To take into account the effect of parton energy loss 
in A A collisions, we use modified effective fragmenta- 
tion functions, Dfi/^{zc,Q'^, L{<f>)) for a produced parton 
c which has to travel in the medium a distance L{(j)) that 
depends on the azimuthal angle in non-central AA colli- 
sions. This will be the dominant source of anisotropy in 
hadron spectra in this parton model. We will use a phe- 
nomenological model 1 20 1 for the modified fragmentation 
functions. The modification in this model depends on 



two parameters: the energy loss per scattering Ec and the 
mean free path Ac for a propagating parton c. The energy 
loss per unit length of distance is then dEc/dx = Cc/Xc- 
We also assume that a gluon's mean free path is half that 
of a quark and then the energy loss dE/dx is twice that 
of a quark. According to a pQCD study of the parton 
energy loss [Ol , 
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where as is the strong coupling constant, Nc = 3 and /i 
is the average transverse momentum kick the propagat- 
ing parton suffers in the medium per scattering. So the 
energy loss per scattering e in our model is now related 
to /z by 
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FIG. 2. Transverse momentum spectra of charged hadrons 
(solid) and tt" (dot-dashed) in central Au + Au collisions 
at ^/s = 200 and 130 GeV without energy loss and with 
dE/dx = 0.5 GeV/fm (dashed lines for charged hadrons and 
dotted lines for vr''.) 

Shown in Fig. g are the calculated spectra for charged 
hadrons, averaged over azimuthal angle, in central Au -\- 
Au collisions at ^ — 130 and 200 GeV with (dashed 
lines) and without (solid lines) parton energy loss. In 
the spectra we also add the soft component in Eq. (H) 
with driah/drj = 520 (700) for ^ = 130 (200) GeV from 
HIJING calculations Q and a^ ~ ^^b'^iax- We used 
a constant dE/dx = 0.5 GeV/fm and A = 1 fm in the 
calculation to demonstrate the effect of parton energy 
loss. The contribution from the soft component is neg- 
ligible at Pt > 2 GeV/c as compared to the hard con- 
tribution when there is no parton energy loss. However, 
when parton energy loss is assumed, the hard contribu- 
tion is strongly suppressed so that the soft component 



become important at some smaller pT- This is why the 
suppression of the large pr spectra due to parton energy 
loss as compared to the spectra without energy loss be- 
comes smaller at smaller p^. However, with dE/dx = 0.5 
GeV/fm, the contribution from the soft component again 
is neghgible at pt > 3 GeV/c. The suppression is pro- 
portional to parton energy loss. The detailed dependence 
of the suppression on the form and values of parton en- 
ergy loss can be found in Ref. [M . Also shown in Fig. 
are the 7r° spectra with (dotted lines) and without (dot- 
dashed lines) parton energy loss. Notice that 7r° spectra 
at high pt are about 40-30% lower than (/i^ + /i+)/2. 
This means that -k"^ only contribute to about 60-70% of 
the total charged hadron spectra. The rest comes from 
kaons and protons. Kaons and protons spectra have dif- 
ferent sensitivity to parton energy loss. This is why the 
suppression factor for large pt tt*' spectra is larger than 
the charged hadron spectra. Detailed study of the flavor 
dependence of the parton energy loss effects can also be 
found in Ref. 11^. 



III. AZIMUTHAL ANISOTROPY 

For non-central ^-|-v4 collisions, the averaged distance 
a parton travels through the medium varies with the 
azimuthal angle and so does the averaged total parton 
energy loss. This will give azimuthal anisotropy in the 
hadron spectra at large transverse momentum with re- 
spect to the reaction plane. The reaction plane is defined 
by the beam direction and the impact parameter b of 
heavy-ion collisions. Experimentally, the reaction plane 
should be determined by low pt particles that consti- 
tute the bulk of matter produced in heavy-ion collisions. 
Large momentum hadrons will not affect the determina- 
tion of this reaction plane. 




r has to travel a distance £(</>, r, b) inside the overlapped 
region in the azimuthal direction <j). With L(0, r, b), one 
can then calculate the modified fragmentation function 
Dhjc{z, Q^ , L{(j), r, b)) and the hadron spectra. The cross 
section is weighted with the overlap function of two collid- 
ing nuclei according to Eq. (|^). To help us to understand 
the azimuthal distribution of hadron spectra, we show in 
Fig. U the averaged distance. 



(i('^)) 



TAB{b) 



d'rU(r)tij(|b-r|)L(0,r,b), (6) 



as a function of </>, that the parton has to travel across 
the overlapped region for different impact parameters, 
where Tab(6) = J d'^rtA{r)tB{\h — r|). As we can see, 
{L{(j))) on the average decreases with impact parame- 
ter. However, the anisotropy increases toward large im- 
pact parameters. Since jet quenching is directly propor- 
tional to {L{(j))), one should expect similar azimuthal dis- 
tribution of the hadron suppression at large pr- The 
azimuthal anisotropy of the hadron spectra, however, 
should depend on both the averaged value of {L{<f))) and 
its anisotropy. 




FIG. 4. The averaged distance a fast parton has to travel 
across the overlapped region of two colliding nuclei as a func- 
tion of the azimuthal angle (j) for different impact parameters. 

The azimuthal angle dependence or azimuthal 
anisotropy of the distance {L{(j))) is directly related to 
the spatial deformation of the dense medium which can 
be characterized by the spatial ellipticity e. For a hard- 
sphere distribution, 



FIG. 3. Illustration of the geometry of the overlapped re- 
gion of two colliding nuclei A and B in the transverse plane. 
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In this model calculation, we will use a hard-sphere nu- 
clear distribution to describe the initial geometry of the 
medium and to estimate the distance of parton propaga- 
tion inside the dense medium. As illustrated by Fig. ^ at 
a given impact parameter 6, a parton produced at point 



is assumed to be given by the initial overlap of two col- 
liding nuclei with a radius Ra ~ 1.12^4^^^ at impact pa- 
rameter 6. Such a simple geometry is far from realistic. 
But it is sufficient to illustrate the qualitative feature of 
azimuthal anisotropy caused by parton energy loss. 



Au+Au (b = 7-8 fm) 




FIG. 5. Azimuthal angle distribution of charged hadrons 
with pT=5,7,and 9 GeV/c in semi-peripheral {b = 7-8 fm) 
Au + Au colhsions at ^/s =130 GeV with parton energy loss 
dE/dx = 0.1(L/fm) GeV/fm. Each distribution is normalized 
to the spectrum in the direction perpendicular to the reaction 
plane. 

Shown in Fig. raare the azimuthal angle distributions of 
charged hadron spectra at different pT normalized to the 
spectra perpendicular to the reaction plane (0 = 7r/2). 
We used dE/dx = 0.1(L/fm) GeV/fm and A = 1 fm in 
the calculation. Due to the azimuthal-angle dependence 
of the averaged total energy loss and the consequent sup- 
pression of large pt hadrons, the hadron spectra show a 
strong dependence on the azimuthal angle. The degree 
of the azimuthal anisotropy is directly proportional to 
the hadron suppression at high pT due to parton energy 
loss. For a fixed value of average total parton energy loss 
the relative effect of parton energy loss on the hadron 
spectra becomes smaller at higher pT |Q. That is why 
the azimuthal anisotropy decreases with pT as shown in 
Fig.|. 

Following a standard procedure in the study of elliptic 
flow [pj , we make a Fourier expansion of the hadron dis- 
tribution in the azimuthal angle. The elliptic anisotropy 
coefficient, V2, is defined as the second order Fourier co- 
efficient, 



V2 



J_ d(j)cos{2(t))da/dyd^pTd(f> 
J_ d(j)dij / dyd^pTdip 
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Shown in Fig. o are the calculated V2 as functions of 
the transverse momentum p^ for different values of the 
parton energy loss dE/dx. The elliptic anisotropy co- 
efficient V2 generally decreases slowly with px but in- 
creases with dE/dx. This is in sharp contrast to the 
hydrodynamic predictions for low pT hadrons as shown 
by the dashed line from Ref. [||. We stopped the cal- 
culation at Pt ~ 3 — 4 GeV/fm for the given dE/dx. 
Below these pT values, the average total parton energy 
loss /S.E = (L) dE/dx will exceed the initial parton trans- 
verse momentum. The hadron spectra from hard pro- 
cesses are strongly suppressed so that contributions from 
non-perturbative processes become dominant and ther- 
mal equilibration and hydrodynamics will determine the 
elliptic anisotropy. Since the coefficient of the differential 



elliptic flow V2{pt) for low px hadrons increases with px, 
the turning point where V2(j>t) starts to saturate and de- 
crease with Pt will provide information on the interplay 
between hard and soft components of particle produc- 
tion. Combined with the magnitude of the azimuthal 
anisotropy at higher px , it also allows one to extract the 
average parton energy loss in the medium. Since one can- 
not make quantitative estimate of parton energy loss, we 
can only give qualitative predictions of V2 and its depen- 
dence on Pt- 
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FIG. 6. The Pt dependence of the coefficient of azimuthal 
anisotropy V2 for charged hadrons for different values of par- 
ton energy loss in semi-peripheral Au + Au collisions at 
\/s — 130 GeV. The dashed line is the hydro result from 



IV. CENTRALITY DEPENDENCE 

In non-central heavy-ion collisions, the average dis- 
tance of parton propagation depends on the impact pa- 
rameter of collisions. The suppression of hadron spectra 
at large pt and the elliptic anisotropy will then also de- 
pend on the centrality. Shown in the upper panel of 
Fig. M is the impact parameter dependence of the calcu- 
lated V2 at fixed pt—5 GeV/c for two different forms of 
parton energy loss, one constant dE/dx = 0.5 GeV/fm 
and another dE/dx = 0.1(L/fm) GeV/fm with a linear 
distance dependence. In both cases, V2 increases with the 
impact parameter, reflecting the increased spatial ellip- 
ticity at large impact parameters. In central collisions, 
the average parton energy loss dE/dx in the second case 
is larger than the constant dE/dx = 0.5 GeV/fm, lead- 
ing to a larger V2. However, as the impact parameter 
increases further, the average parton energy loss per unit 
length with a linear distance dependence becomes smaller 
than the constant dE/dx = 0.5 GeV/fm. This leads to 
smaller V2. Eventually, in more peripheral collisions, V2 
decreases with the impact parameter in the second case. 
Therefore, as compared to the case of a constant dE/dx, 
a linear distance dependence of dE/dx gives a slower im- 
pact parameter dependence of V2 . 

Since the elliptic anisotropy in hadron spectra is di- 
rectly related to the spatial ellipticity, it is useful to study 



the impact parameter dependence of the ratio f2/e as 
shown in the lower panel of Fig. [7|. Such a ratio will 
also reduce the sensitivity to the modeling of the ini- 
tial density distribution of the overlapped region, e.g., 
hard sphere vs. Woods-Saxon nuclear distributions. For 
both forms of parton energy loss, the normalized ellip- 
tic anisotropy coefficient increases with impact parame- 
ter very rapidly in central collisions and then starts to 
decrease at relatively small impact parameters. The de- 
crease is much faster for the case of linear distance depen- 
dence of dE/dx. This feature is markedly different from 
the hydrodynamic models |p[, in which V2/e remains con- 
stant for a very large range of impact parameters. So far 
experimental measurements of the impact parameter de- 
pendence of V2 [^f y for low pt hadrons show a behavior 
roughly similar to the hydrodynamic results. A different 
behavior of V2/S at large pr that decreases with the im- 
pact parameter will be a clear indication of the effect of 
parton energy loss. 
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FIG. 7. Upper panel: The impact parameter dependence 
of azimuthal anisotropy coefficient V2 for charged hadrons 
with pt~5 GeV/c in Au + Au collisions at y^ = 130 GeV 
for a constant dE/dx = 0.5 GeV/fm (open circle) and dis- 
tance-dependent dE/dx — 0.1(I//fm) GeV/fm (solid circle) 
parton energy loss. Lower panel: The impact parameter de- 
pendence of V2/e. Here, e = 6/27?a is the spatial ellipticity of 
the initial overlapped region of collisions with a hard sphere 
nuclear distribution. 

For a qualitative study in this paper, we will neglect 
the effect of expansion and we are only concerned with 
the average parton energy loss. In this case, the parton 
energy loss dE/dx should depend on the averaged par- 
ticle density in the overlapped region of heavy-ion col- 
lisions, which in turn depends on the centrality of the 
collisions. 

According to Eq. (0), dE/dx is proportional to /i^/A — 
/i^crp with /i being the momentum transfer per scat- 



tering. For a parton scattering cross section a that is 
screened by /j,^, /i^tr is roughly a constant. Then dE/dx 
will be proportional to the parton density p. We as- 
sume that the initial parton density is proportional to 
the final hadron multiplicity per unit transverse area, 
p ^ dN/dy/{ToS{b)), with S{b) the transverse area of 
the initial overlapped region and tq the initial time. Us- 
ing the two-component (soft and hard) model of particle 
production in HIJING |l9], dN/dy has one term (soft) 
proportional to the number of participant nucleons A^part 
and another (hard) proportional to the number of bi- 
nary collisions A^binary In a Glauber model of nuclear 

coUisions, S{h) ^ Nlilih) and A^binary(6) ~ Kilih). 
Parameterizing the centrality dependence of dN/dy in 
the HIJING calculation |gj] for Au + Au collisions at 
-yi = 130 GeV, we have 



dE 
dx 



rl/3 



pi{b) - 0.072N^i:,ib) 1 + 0.13N^i:,ib) 



ri/3. 



(9) 



where the impact parameter dependence of Npart{b) can 
be calculated in a Glauber model of nuclear collisions 
with Woods-Saxon nuclear distributions. The above 
centrality dependence of the average dE/dx is normal- 
ized to central collisions, dE/dx{b = 0) = {dE/dx)o, 
or p(b = 0) = 1. Such an impact parameter depen- 
dence of dE/dx is plotted in Fig. g as the solid line for 
{dE/dx)o = 0.1(L/fm) GcV/fm. 
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FIG. 8. Three scenarios of impact parameter dependence 
of parton energy loss dE/dx — 0.1{L/im)p{b) GeV/fm 

Theoretical studies have predicted a large parton en- 
ergy loss in a hot partonic medium as compared to cold 
nuclear matter ||l^. However, analysis of the large pr 
7r° spectra in central Pb + Pb collisions from WA98 ex- 
periment Ig^ at the SPS energy V4 = 17 GeV shows no 
indication of parton energy loss ||26| . Since the high pT 
spectra at SPS are very sensitive to any change of the 
jet cross section, the constraint on parton energy loss by 
the WA98 data is very stringent and gives a limit that 
is much smaller than the most conservative estimate of 
parton energy loss in a hadronic matter [ 16| . This casts 
serious doubts on the accuracy of the theoretical esti- 



mates which are based on a scenario of a static and in- 
finitely large (L/A ^ 1) dense partonic gas. Though 
recent studies p7|-p9| have considered the finite number 
of scatterings in a medium with finite size, it is still not 
clear whether the estimates are numerically consistent 
with the WA98 data. The experimental data of WA98 
could also imply that the lifetime of the dense partonic 
system is too short to induce any parton energy loss. 
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FIG. 9. Impact parameter dependence of V2 and v^je for 
charged hadrons with pt=5 GeiV /c in Au + Au collisions at 
y/s = 130 GeV for three different scenarios of density depen- 
dence of parton energy loss. 

At RHIC energies the initial energy density is expected 
to be higher than at the SPS with longer lifetime. Indeed 
the first experimental measurement of dNch/drj by PHO- 
BOS 1^21 shows that the particle production for central 
Au + Au collision at ^/s = 130 GeV is about 50% higher 
than in central Pb + Pb collisions at the SPS energy, 
consistent with the default value of HIJING simulations. 
If we take the implication of the WA98 data at its face 
value, we can assume that the effective parton energy 
loss vanishes in non-central collisions in which the aver- 
age particle density is equal to or less than that of cen- 
tral Pb + Pb colhsions at the SPS energy. Using HIJING 
simulations we find that the particle rapidity density per 
unit transverse area in \ri\ < 1 for Au + Au collisions be- 
comes less than that of central Pb + Pb collisions at the 
SPS energy for impact parameters larger than 7 fm at 
Vs = 130 GeV, or 8 fm at Vs = 200 GeV. For smaller 
impact parameters, we will assume the normal centrality 
dependence of the parton energy loss in Eq. M) . Such a 
scenario is denoted as density dependence pii{b) and is 
shown in Fig. @ as the dashed line. The sudden drop of 
dE/dx with impact parameter might be extreme and is 
only to illustrate the effect of a possible critical behavior 
in jet quenching. 

Shown in Fig. are the calculated coefficients of the 



azimuthal anisotropy V2 and V2/£ for charged hadrons 
at pt=5 GeV/c with the above two scenarios of den- 
sity dependence of parton energy loss (solid and dashed 
lines). As a comparison, we also show the results of a 
density- independent dE/dx (dot-dashed lines). Compar- 
ing the results of scenario I of density dependence (solid 
lines) to a density-independent dE/dx (dot-dashed lines), 
it is obvious that the density dependence of the parton 
energy loss dE/dx reduces the azimuthal anisotropy at 
large impact parameters because of the reduced parton 
energy loss in peripheral collisions as compared to the 
central ones. The ratio V2/S therefore decreases faster 
toward the peripheral collisions as compared to the case 
of a density-independent dE/dx. Because of the com- 
plication of such a density dependence, it is difficult to 
single out the quadratic distance dependence of the par- 
ton energy loss by studying the centrality dependence of 
the azimuthal anisotropy alone. One has to combine it 
with the study of density dependence by measuring the 
anisotropy for a fixed system at different colliding ener- 
gies. In the second scenario of the density dependence of 
parton energy loss, V2 will vanish above impact parame- 
ter 6 > 7 fm, where dE/dx is assumed to vanish when the 
particle density is below what was reached in the central 
Pb + Pb collisions at the SPS energy. The drop should 
occur at larger impact parameters for higher colliding 
energies because of the increased particle density in the 
collision region. The behavior of V2 at small or interme- 
diate impact parameters is dictated by the geometry of 
the medium and is not very sensitive to small changes of 
the density dependence of the parton energy loss. It is 
more sensitive to the density dependence at large impact 
parameters in more peripheral collisions. The vanishing 
dE/dx at large impact parameter in the second scenario 
has very dramatic effect on V2- Such a scenario will be 
easy to verify from experimental data. 

Since the azimuthal anisotropy in large px hadron 
spectra in a parton model is caused predominantly by 
jet quenching or parton energy loss, such a study should 
be complemented with direct measurements of the sup- 
pression of large px hadron spectra relative to pp or pA 
collisions at the same energy. To facilitate such direct 
measurements, one needs to scale the spectra in A A col- 
lisions by that of pp and geometric factors. According to 
the parton model in Eq. (fil), the ratio 



RAAm,PT) 



dNAA/dyd'^PT 



{Nbin^,y){{b))dNpp/dyd^PT 



(10) 



will be 1, independent of px and impact parameter b if 
there are no nuclear effects in the large pr hadron pro- 
duction. Jet quenching will suppress the large pr spectra 
and reduce the above ratio for non-vanishing parton en- 
ergy loss Q. Here (A'^binary )((&)), defined as 



(A^binary)((fo)) = 



_PP 



^AA 



[b 



min-i ^maxj 



d%TAA{b), 



(11) 



is the averaged number of binary collisions that can be 
calculated in a Glauber model of nuclear collisions fig ] 



with given (b) ~ [bn 



.)/2. 
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FIG. 10. Impact parameter dependence of the ratio 
Raa(pt, h) of the charged hadron spectra at pt=5 GeV/c in 
Au + Au and pp collisions at y^ = 130 GeV for three different 
scenarios of density dependence of parton energy loss. 

Shown in Fig. |l^ are the calculated ratios for charged 
hadron spectra at pr = 5 GeV/c in Au + Au collisions 
at y/s = 130 GeV as functions of the impact parameter 
with three different density dependencies of parton en- 
ergy loss. The suppression due to parton energy loss in 
general will decrease with impact parameter since the size 
of the overlapped region decreases toward more periph- 
eral collisions. For very peripheral collisions, one should 
recover the results of pp collisions and the ratio should 
be 1. However, if the parton energy loss dE/dx depends 
on the average particle density, the ratio (solid) should 
approach 1 more rapidly than in the case of a density- 
independent dE/dx (dot-dashed). Again, if we assume 
that parton energy loss vanishes at 6 > 7 fm from con- 
straints by WA98 data at the SPS energy, the ratio will 
become or is larger than 1 already in semi-peripheral col- 
lisions. The ratio (dashed) will actually become larger 
than 1 because of initial px broadening or Cronin effect. 
Such a centrality dependence of hadron suppression at 
large pt combined with the energy dependence will shed 
light on parton energy loss in the medium. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper we have studied in a parton model the 
azimuthal anisotropy of large pt hadron spectra in non- 
central high-energy heavy-ion collisions due to parton en- 
ergy loss. We demonstrated that the anisotropy is very 
sensitive to the parton energy loss dE/dx and its depen- 
dence on the size and density of the medium. We predict 
that the coefficient of the anisotropy V2 decreases slowly 
with Pt- The parton model also predicts an early satura- 
tion of V2 as a function of the impact parameter b. Once 
divided by the spatial anisotropy e(6), the ratio V2/ s will 
increase with b for very central collisions and then will 



decrease toward peripheral collisions. This is in sharp 
contrast to the centrality dependence of the elliptic flow 
V2 for low Pt hadrons. Hydrodynamic models & predict 
almost a constant V2/ s in a very large range of central- 
ities. We found that the centrality dependence of V2 is 
sensitive to the size and density dependence of parton 
energy loss dE/dx in the medium. 

Since the azimuthal anisotropy is directly related to 
parton energy loss, one should combine the study of 
azimuthal anisotropy with the direct measurements of 
hadron spectra suppression at large pt- We propose to 
study the density dependence of parton energy loss by 
measuring the hadron suppression and the anisotropy at 
different colliding energies for the same system. The 
centrality dependence of the hadron suppression and 
anisotropy can then provide information on the size de- 
pendence of the parton energy loss. Perturbative QCD 
studies predict a nonlinear dependence of the total en- 
ergy loss due to non-Abelian gluon radiation. 

Based on constraints by WA98 experimental data on 
high Pt pion spectra, we also proposed a density depen- 
dence of the parton energy loss that will vanish in pe- 
ripheral collisions, where the particle density is equal to 
or smaller than what has been achieved in the central 
Pb + Pb collisions at the SPS energy. Then both hadron 
spectra suppression and the azimuthal anisotropy will 
disappear in these peripheral collisions. Such a strong 
centrality dependence would indicate an onset of parton 
energy loss in a dense medium. 

Parton energy loss is always associated with momen- 
tum broadening perpendicular to the parton propaga- 
tion direction. Such transverse momentum broadening 
can also lead to azimuthal anisotropy in the final hadron 
spectra. One can characterize the broadening by the mo- 
mentum transfer /.t in each parton-medium scattering. 
The effect of the broadening on the azimuthal anisotropy 
can then be characterized by the diffraction angle of each 
scattering 0o = /i/pT- For dE/dx = 0.1(L/fm) GeV/fm 
and A ~ 1 fm, which we have used to demonstrate the 
effect of parton energy loss in this paper, /i ~ 0.3 GeV/c 
according to Eq. (||) and (/)o ~ 0.06 for p^ — 5 GeV/c. In 
a medium with a size of only a few mean free paths, the 
contribution to V2 due to such a small momentum trans- 
fer is almost negligible. Therefore, the dominant cause of 
elliptic anisotropy at large px will be the radiative parton 
energy loss in the medium. If any suppression of large 
Pt hadron is observed in experiments, there should be 
non-vanishing azimuthal anisotropy, and vice versa. 

We have considered a parton energy loss that has only 
a weak (logarithmic) dependence on the initial parton 
energy. Recent studies p7| , however, show a stronger 
energy dependence as a result of the kinetic limits in the 
estimate of energy loss for a parton with finite initial 
energy. Such a strong energy dependence would give a 
different px dependence of hadron spectra suppression 
and the accompanying azimuthal anisotropy. 

For a qualitative study in this paper, we did not con- 
sider the expansion of the system. The longitudinal ex- 



pansion will change the particle density in the medium, 
while the transverse expansion will decrease the spatial 
anisotropy as the system evolves with time. The decrease 
of particle density will reduce the average parton energy 
loss over the course of the evolution. The change of spa- 
tial anisotropy will reduce the resultant anisotropy in 
hadron spectra. It could also affect the centrality de- 
pendence of the jet quenching effect. Detailed study of 
the jet quenching in a dynamically evolving system will 
be the subject of future studies and is needed for a more 
quantitative analysis of any experimental data. 
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